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Abstract 

We show that there are inflationary models for which perturbations in the energy 
momentum tensor, which are of second order in the scalar field, cannot be neglected. 
We first specify the conditions under which the usual first order perturbations are 
absent. We then analyze classically, the growth and decay of our new type of pertur- 
bations for one mode of fluctuations Stfik in the scalar field. We generalize this analysis, 
considering the contribution from the whole spectrum of 5(f) to a given wavelength of ge- 
ometrical perturbations. Finally, we discuss the evolution of the perturbations during 
the subsequent radiation dominated era and discuss the resulting spectrum of density 
fluctuations. In the case of a massless scalar field we find a spectral index n — 4. For 
massive scalar fields we obtain n — but the resulting amplitude of fluctuations for 
inflation around a GUT scale are by far too high. Hence, 'conventional' inflationary 
models must not be influenced by this new type of perturbations, in order to lead to 
acceptable perturbations. 



PACS: 98.80.Cq 98.65.D 



1 Introduction 



Inflation was originally proposed |lj |2| to solve some shortcomings of the standard cos- 
mological model. In addition, inflation can provide initial perturbations needed for the 
formation of the structures observed in our universe. Most inflationary scenarios predict a 
scale-invariant spectrum of fluctuations §, known as Harrison-Zel'dovich spectrum |||, ||], 
which was found to be consistent with observations performed by the COBE satellite |], 0]. 
This remarkable agreement motivated many of us in studying the evolution of perturba- 
tions in the context of inflation. In particular, since the amplitude of fluctuations poses 
the stringest constraints on inflationary models , we are interested in analyzing the growth 
and decay of cosmological perturbations. 

To our knowledge, up to now the only perturbations of the energy momentum tensor 
which have been considered, were first order in 8(f). (Recently we where told that quantum 
fluctuations which are second order in 6(f> have been investigated ||.) On the other hand, 
since successful inflation requires that 4> and V,^ are small enough, one can envisage situ- 
ations where <p is of the order of 5(f) and V,<p <C Vm 8(f)- In this case, all contributions to 
the perturbation of T^ u are second order in 5(f). However, it has to be stressed that these 
perturbations correspond to first order perturbations in both, T^ u and the geometry, and 
thus they can be analyzed by usual cosmological perturbation theory. This is the issue 
addressed in this paper. 

Our work is organized as follows: In the next section, we analyze the conditions on 
the potential for which first order perturbations can be neglected. In section 3 we study 
the growth and decay of geometrical fluctuations induced by a single mode perturbation 
in 5(j). We proceed in section 4 by generalizing to a generic spectrum of 8(j). In this case 
the whole spectrum contributes to the geometrical perturbations of a given wavelength. 
In section 5, the evolution of perturbations is followed through the radiation dominated 
era until second horizon crossing to calculate the resulting power spectrum of density 
fluctuations. We finally summarize our conclusions. 

Since we are only concerned with the very early universe, we neglect spatial curvature. 
Furthermore, we use the signature (— , +, +, +), such that in conformal time rj, the metric is 
given by ds 2 = a 2 (—dr/ 2 + dx 2 ). We denote by prime and dot the derivatives with respect 
to conformal and cosmological time respectively. In what follows, G is the Newtonian 
gravitational constant, mpi = 1/VG is the Planck mass (c = U = 1) and H = d/a = a' /a 2 
is the Hubble parameter. 



2 When are second order perturbations important? 

Let us consider a minimally coupled real scalar field (f) with Lagrangian density 

L= X -d^4>-V{4>). (l) 

The stress energy tensor of (j) is given by 

T^v = d^(j)d v (j) - Lg^ u . (2) 

During inflation, the energy density of the universe is dominated by the potential term V 
and thus assuming a homogeneous and isotropic inflaton field (J), we have (l/2)(p 2 <C V((j)). 
The equation of motion for <p(t) 

4> + 3H0 + ((/)) = (3) 
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(where V,* denotes the derivative of V with respect to <f>) and Friedmann's equation, 
H 2 = (8vrG/3)F, imply |V, | < V48nV/m P i. 

First order perturbations in the energy momentum tensor @ have the form 

j>5<j> or V )4 ,5<f>. (4) 

We now want to study under which conditions the first order perturbations in 8(f> are 
negligible, namely, when <f) ~ 6(f) and 

\V^\<x^{\V,^S<l>\,\-5(f>%\\5f\,^\A{Sif>)\} , (5) 

Ob Ob (X 

where A denotes the usual three dimensional Laplacian. Let us consider the scalar field 
to be of the form 

4> = 4>o + Sep, with 4>q = constant and 5(f) = 5(pk = fnA- k {rf)e lkx ' , (6) 

where we cast the dimension of 6(f) in the mass parameter m, so that A k is a dimensionless 
quantity. In the next section we will appropriately choose the value of m. Provided 
condition (]5|) is met, we can then write the equation of motion for 5<p as 

Al + 2HaA' k + A k (k 2 + a 2 B) = 0, (7) 

where 

B = V, H (0o )• 

We want to emphasize that the geometrical terms usually present in the above equation 
(see, e.g. ||]) do not contribute in our case, up to first order, since and are already 
first order perturbations. The solution of (|7|) is given by 



A k = A in (k)y^ 2 J±„(y), with V = ^ and ^^y 1 "^' ( 8 ) 

where by J a we denote the Bessel function of order a. As we will show later, for suitable 
initial conditions Ai n {k) is actually independent of k, so that it can be absorbed in the 
definition of m. Thus, without loss of generality we set Ai n {k) = 1. For the space and 
time derivatives of 5(p we then find 

W = -^W /2 [3J±v(y)+yJ±u-i(y)-yJ±u+i(y)]e ikx , (9) 

S<f> }j = ik jm y^ 2 J ±u (y)e ikx . (10) 

Before proceeding with the discussion of the conditions imposed on the form of the po- 
tential V, we will distinguish between super- and sub-horizon perturbations. (Here, as 
usual, super (sub)- horizon denote larger (smaller) than the Hubble radius respectively.) 

For super-horizon perturbations, A = a/k S> 1/H, i.e., y <C 1. For the least decaying 
mode we then find 

\6<j>\=H\(l-v)8<j>\. (11) 
For B > 0, 

|^ 2 | < B\8(P\ 2 « V&o), 
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where the last inequality sign has to hold in order for the universe to be in an inflationary 
phase. Therefore, if \5(p\ was small initially, it remains small during inflation for arbitrary 
positive values of B (which is of course also intuitively clear). Furthermore, on super- 
horizon scales |A5</>| = /c 2 |<50| <C \5<t>"\ and the last term in (5) can thus be neglected. The 
limit on V,^ from (|^) then reduces to \V,^ ((fto)\ ^ -B|50|. After iV e-foldings of inflation 
\5<fi\ ~ m\ exp[N(u — 3/2)]|, so finally the conditions imposed on V are 

|V,0(0o)|< 5mexp(-3iV/2), for B > H 2 (12) 
|V, (0o)|< Bmexp[-N(%-i>)], for B < H 2 (13) 

On the other hand, for B < 0, grows according to oc a u ~ 3 ^ 2 . For 5cj) not to grow 
extensively during N e-foldings of inflation, we hence have to require 

N(y - 3/2) < 1 or, equivalently, - H 2 /N < B. 

In that case the conditions on V become 

— H 2 

\V,f (0 O )| « Bm and B > — . (14) 
For sub-horizon perturbations, A = a/k -C 1/H, i.e., i/ > 1, we have 

1^1 w mHy 2 w Hy\S<f>\ = -\5S (15) 

a 

and 

\8<f>,j | « kA8<j)\ ps a\6j>\. (16) 
Therefore, the conditions on the potential V result in 

|V, (0 O )| «max{ J B|^|,^|^|}. (17) 

On sub-horizon scales, classical perturbations decay and quantum effects are important, 
and a quantum analysis has to be performed. 

Provided < B <C H 2 , the < (5<t>) 2 > can be approximated by 

which suggests that [[R]] 

_ a 3 gVB + (W fc 2B/3g2 _ y/l + B/(yH)* 2B/3H 2 
Uk ~ 2fc3 W 2y2 V 

denotes the number density of c/>-particles. On sub-horizon scales, <C 1, which means 
that the quantum analysis is needed. On the other hand, on super-horizon scales, n& S> 1 
and thus our classical analysis is indeed sufficient. 

A quantum analysis study in the case of sub-horizon perturbations shows that the 
main contribution to (|i~8| ) is given by the second term, which approximately yields 
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A natural choice for m, in the case that B <C H , is therefore m = H/(2ir), leading to 

x)= r ^(4-)^ju4z)d% 







2vr y Ha' y Ha' 



for super-horizon perturbations. 

In the situation where the above imposed conditions on the potential V, namely equa- 
tions (|l2|), (p^|), and ([l?]) are satisfied, the perturbations in Ij^ are given by 

ST^ = (<ty),„ (6</>), v -\g^{5(f>),x {5(f)), X +9^V, H 6<f>6(j). (22) 

In the following two sections, we will study the behavior of the geometrical pertur- 
bations induced by the above 5T 11V . Before proceeding with this study, we would like to 
briefly discuss some inflationary models, for which our conditions are met. 

• Old Inflation: In old inflation |l]], the inflaton field (f> is sitting in a minimum of 
the potential V, so that B > and V,s = 0, which means that our conditions are 
satisfied. However, since in this model the phase transition is of first order and, 
thus, proceeds by bubble nucleation, it has the problem of the so-called graceful 
exit |H| , |i"2"| . This means that the bubble nucleation process is slow with respect to 
the universal expansion and, therefore, the phase transition can never be completed. 
On the other hand, if this problem is solved, (e.g., by extended or hyper-extended 
inflation fl3| ) then the density fluctuations produced by bubble nucleation will, most 
probably, be the dominant ones. 

• New Inflation: In new inflation [jEj], the inflaton field (f> is sitting either in a 
minimum or in a maximum of the potential V, depending on the details of the 
model. Nevertheless, if B < 0, the slow rollover condition requires \B\ to be small. 
Thus, our condition, namely, B > —H 2 /N, may well be satisfied. In such a case, one 
might think that the problems usually arising with first order perturbations (i.e., 
the requirement of extremely flat potential and of a very weakly coupled scalar field) 
can be avoided. However, a more detailed study has to be performed. 

• Supersymmetric Inflation: In certain models of supersymmetric inflation [jljj], 
both V,0 and B vanish and our conditions are certainly satisfied. Thus our analysis 
can indeed be applied. 



Chaotic Inflation: In chaotic inflation [16|, V oc \<fi n , so that V,<f, ~ B(f>o and 
therefore, V,^ ^ BS(p requires that (fio ^ 5(f), which is of course inconsistent with 
perturbation analysis. As a conclusion, within the context of chaotic inflation our 
second order perturbations are irrelevant. 

de Sitter universe: Regarding an arbitrary scalar field within an inflationary 
model, our second order perturbation analysis can be applied to study the fluc- 
tuations of that field, provided either the potential is very flat, or, the field is sitting 
in a minimum of a potential, which then can be of an arbitrary form. 



3 Perturbation analysis for one mode 

In this section we will investigate the evolution of second order perturbations in the scalar 
field, which are of the specific form 5(f) = mAk(rj)e l]i ' x . Since the quantum analysis is only 
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valid for B <C H 2 , we now restrict ourselves to only this case. As we have seen before, 
the classical analysis applies only for super-horizon perturbations and the value of 5(f) at 
horizon crossing is given by 5<p{y = 1) = H/(2tv). The solution, Eq. (^), of the classical 
equation of motion, Eq. (|7|), for 5(j) implies that the matching condition at horizon crossing, 
with m = H/(2n), yields 

A in (k) = 1/J±„(1) ~ 1 , 

which is /c-independent. This justifies to drop the prefactor Ai n {k) in the sequel. We now 



5p = i6Ti = ^f p , (25) 



parametrize the energy momentum tensor Eq. (22) by 

5T 00 = m 2 k 2 f p , (23) 
ST 0j = m 2 kf v , 1 , (24) 

3 1 a 2 

S Tj i = m 2 k 2 (f 7T , j i-^Af 7T 5 j i), (26) 

where 8rji is the traceless part of the perturbed stress tensor. In terms of the function A k 
defined above, we obtain 

f P = ^[A' 2 -A 2 (k 2 -a 2 B)]e 2 ^, (27) 
f P = ^[\k 2 A 2 + A' 2 -a 2 BA 2 ]e 2 ^^ (28) 

fv = ^A k A' k e 2 ^, (29) 

U = A 2 k e 2ik *. (30) 

Using the perturbed Einstein equations, we will now relate this energy momentum ten- 
sor to the perturbations in the geometry. In the case of scalar perturbations, which are the 
relevant ones for structure formation, geometrical perturbations can be expressed in terms 



of the Bardeen potentials $ [17]. In this paper we only analyze scalar perturbations. 

Let us now briefly discuss the physical meaning of the Bardeen potentials. They were 
originally introduced because they are invariant under linearized coordinate transforma- 
tions, which, in the context of cosmological perturbation theory, are called gauge trans- 



formations. As it has been shown [18, [T^] (Stewart lemma), gauge invariant quantities are 
perturbations of tensor fields with either vanishing, or constant, background contribution. 
Since homogeneous and isotropic universes (Friedmann-Lemaitre models) are conformally 
flat, they have a vanishing Weyl tensor. Thus, according to the Stewart lemma, its per- 
turbations are gauge invariant. The Bardeen potentials are related to the electrical part 
of the Weyl tensor by 

Eij = a- 2 C l0j0 = \[($- ~ \H® ~ *)] ; (31) 

the magnetic part vanishes for scalar perturbations pCfl . The amplitude of perturbations 
in the geometry is typically given by 

A- 2 

,K ',/\/7/r;n !■ ~ - ' 
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where Rg is the Ricci curvature. In the inflationary universe 
A^y 2 {<$>-^) . 



(32) 



Hence, the Bardeen potentials determine the amplitude of perturbations at horizon cross- 
ing. On super-horizon scales, where our analysis is relevant, the amplitude of perturbations 
is constant only if the Bardeen potentials grow as a 2 . 

Within the context of this paper, 5T^ represents the perturbation of the constant tensor 
field y(0o)<^, thus it is gauge invariant and it can be regarded as a seed perturbation 
according to The Bardeen potentials are then given by [21] 



-e[f p k 2 + 3(a'/a)kf v ] ; 
-2eAf n , 



(33) 
(34) 



where e = A-K(m/mpi) 2 must be much smaller than 1 for linear perturbation analysis to 
hold. For B -C H 2 , e becomes e = H 2 /(irm'p l ). Equations (|33| ) and (34) imply that in our 
case the Bardeen potentials are 



e 



4? 



A\(k 2 -a 2 B) + 3-A k A[ 



2ikx 



e 

2F 



A' 2 + A 2 (3k 2 + a 2 B) + 3-A k A' k 

a 



„2ik-x 



On super-horizon scales (y < 1), A k given by Eq. (||) reduces to 
2 U 



A k 



T(l-u 



-V 



3/2-i> 



Inserting this result in equations (£3^) and (36), we find 

_ 2 2u-l 



[r(i - u)] 2 



3-2v 2ik-x . 



o . o2i/-l 

z e „,3-2i/„2ikx 



-y 



[r(i - u)f 

For the relevant regime, B <C H 2 and the above expressions lead to 



3$ = * 



-3 e 



7T 



,2ikx 



(35) 
(36) 

(37) 

(38) 
(39) 

(40) 



The Bardeen potentials are constant in time meaning that the amplitude of curvature 
perturbations decays exponentially. 

We shall see in Section 5, how the amplitude e of \& is connected with observed fluc- 
tuations in the cosmic microwave background and what its value has to be due to a 
normalization, e.g. according to the COBE quadrupole. 



4 Perturbation analysis for the whole spectrum 

So far we have only considered the contribution of one wave vector k. However, since the 
Bardeen potentials are quadratic in 5(p, to a given wave vector ki of the (<!>, ^fr) spectrum, 
there is a contribution from all wave vectors k of 5(p. We therefore address in this section, 
the contribution of the whole spectrum of k in the evolution of the Bardeen potentials. 
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We use the general ansatz 
5<j) = mA(x, if). 

The fp,f v ,f n functions, which enter in the Bardeen potentials, become 

A(x) = A 2 ; (41) 

Vjf v (x) = lA'VjA- (42) 

/„(x) = l 2 [A' 2 + (VA) 2 + a 2 BA 2 ]/2 , (43) 

where I denotes an arbitrary length scale, for example I = l/k±. To find the spectrum of 
the above functions we now perform a Fourier transform. To obtain a finite result, we have 
to introduce a cutoff. Taking into account that the fluctuations A(k) are distributed with 
uncorrelated phases, it makes sense to neglect contributions from wavelengths A < Ai, 
i.e., to choose the cutoff k\ and correspondingly the normalization volume kf. Using the 
convolution theorem, we then obtain 

1 r kl ,s ; 



A(kl) w J^]wJo d6kA (\M)A(\^--k\); (44) 
/«(ki) ~ (27r) 3 /2fc 4 <PkA'(\k\)A{\ti - k|)kj • (k! - k) ; (45) 

/pCkx) - ^^r^'(|klK(|k 1 -k|) + 



(k • (kx - k) + a 2 B)A(\k\)A(\k 1 - k|)] . (46) 



With 



• 4 (' k " = fo^)^'"' 

the above integrals finally yield 



M kl) - Vg^ ; (47) 

/p(ki) - f^ + hi*]*^, ( 49 ) 



where yi = k\/{Ha). Inserting these results in equations (|33| and 34), we find for the 
Bardeen potentials 



d> „ ~ e ^ . (50) 
5 

This result shows that the Bardeen potentials remain constant during the inflationary 
period and the geometrical perturbation amplitude A, introduced in the last section, thus 
decays like 1/a 2 . However, the relevant physical quantity for structure formation, is the 
change of 5p/p between the first and second horizon crossing. At horizon crossing, 5p/p is 
well approximated by the Bardeen potentials 0. Hence, we will now study the evolution 
of the Bardeen potentials during the radiation dominated era, until the time of second 
horizon crossing. 
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5 The radiation dominated era and second horizon crossing 



In this section H without argument denotes the value of the Hubble parameter during 
inflation. The value of a quantity at the end of inflation is indicated by the subscript end 
and at present time by a subscript o- 



5.1 The Bardeen potentials at second horizon crossing 

According to the general equation of motion for S(f>, (j^), the evolution of the perturbation 
amplitude in the radiation dominated regime is given by 

^ + -^ + (lWM fc = 0, (51) 
dy z y dy 

where y = k/(aH(y)) = krj in the radiation dominated epoch, and a = a 2 B /(k 2 rj) 2 = 
const. Since at the end of inflation rj = X/{aH) to a very good accuracy, we find 

B 

a 



Here y end is the value of y at the end of inflation and H denotes the Hubble parameter 
during the inflationary period. For B = 0, ( |5ll) is exactly solved in terms of Bessel functions 
and the 'growing mode' is given by 

A k (y) = —Ji/ 2 (y) ■ 

This solution is still approximately valid as long as ay 2 < 1. To require oty 2 nd < 1 for all 
physically interesting scales, Hao < k < HqUq yields in terms of the mass parameter B 

B « z 2 end H 2 ~ (lCT 6 ey) 2 • (T end /W 15 GeV) . (52) 

Obsevationally, such a small mass can not be distinguished from B = 0. The general 
solution to (51) at y = 1 can be found by exponentiating the matrix 

M(y,y end ) = ( _ {1 l ay2) _l /y )dy 



yend 



( y-yend \ 

\-(y- Vend + a/3(y 3 - yl nd )) 2\n(y end /y) ) ' 

which determines the solution according to 

(S)^ p(M( " J, ("fc))- 

Unfortunately M cannot be diagonalized; but explicit exponentiation (by expansion in 
power series) yields a simple, summable result if \M2i\ 3> |-/Wi2|> IM22I, which is valid for 
a S> 1 and y 3> y en d, In the limit a>lwe find 



exp(M(l,y end )) w (1 + ln(y end )) 



cos y/a/3 y/3/a sin \/a/3 
- \/a/3 sin \/a/3 cos \/a/3 
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For each element in exp(M) we have kept only the highest power in a. Using the exact 
solution during the inflationary epoch, eqn. (37), which yields 

A'(y end ) = ^yl nd A(y end ) , (53) 

we obtain for A k , A' k at y = 1 

= (1 + ln(y en d))[cos \f*/3 + yl nd \[a/?> sm. \[a/3]A k (y end ) 

A'ki 1 ) = 0- + ln (2/end))[(a/%end cos \f®/3 ~ \fa/3sm \fa/3]A k (y end ) . 
In the radiation dominated era eqn. (|36| ) for the Bardeen potential \& yields 

= -(^ 2 )[(^) 2 + ^(3 + ay 2 ) + (3/l/Mp . (54) 
For q»1 (which implies y enc j <C 1) this yields 

*rad(l/ = 1) ~ - e -^A 2 {y end ) . (55) 
Equation (||) at y = y end with the initial condition (|53|) leads for ay 2 nd S> 1 to 

^radiVend) ~ -eay;? rad ,4 2 (y end ) tx y~ n 2 d . 

On the other hand, from the inflationary universe calculation we have ^infiUend) = e = 
const. Assuming now, that the Bardeen potentials do not jump when the universe enters 
the radiation dominated phase, \P has changed since the end of inflation until second 
horizon crossing by a factor 

9(k) = ^rad{t)/^rad{yend) ~ Vend/ 3 • ( 56 ) 

This yields for the final Bardeen potential at second horizon crossing 

*(Vk)*-^ = - <H ^ d)2 for a = B/(HV end ) = BH 2 ai nd /k 4 »l, (57) 

where r\ k denotes the time of second horizon crossing of wave number k, i.e. rjk = 1/k and 
e = H 2 1 '(irrripi) . In the case B = 0we find ^ ra d ~ const., independent of y end . This leads 
to a growth factor of g(k) ~ 1, and thus 

~ -e for 5 = . (58) 



In the intermediate regime, O(a) ~ 1, we have solved ( pl|) numerically. The solutions A k 
for some values of a are represented in Fig. 1. 

In Fig. 2 we plot the value of ^f(rjk) m units of e as a function of physical wave number 
k/a in units of H (In other words, we plot f asa function of y en d-)- We chose the parameter 
value B/H 2 = 0.01. For small wave numbers, q>1, the amplitude of ^ decays like k 
until approximately ay 2 nd = B H 2 / (k / a end ) 2 ~ 1. At (k/a end ) 4 ~ BH 2 , i.e. a ~ 1 the 
potential bends over to a constant. 

This results remain qualitatively valid also in the matter dominated era which we thus 
shall not treat seperately. 
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5.2 The power spectrum 

Let us now briefly clarify what power spectrum results from the findings of the last subsec- 
tion. Since the explanations of spectral index and scale invariance found in the literature 
are often somewhat confusing, we shall be rather explicit in this paragraph. 
The power spectrum is defined by 

P(k, v ) = \5 k (i 1 )\ 2 , 

where S k (r]) is the Fourier transform of the density fluctuation (8p/ p)(x,rj). Clearly for 
perturbations which are larger than the size of the horizon, k < 1/rj, P(k,r]) is gauge 
dependent. But we are mainly interested in P(k, t]q) on scales k > l/r/o and a^/k ~ Ao > 
lOMpc, i.e., scales which are small enough to be observable today and large enough, so that 
they are probably not severely affected by nonlinear clustering. The spectrum P(k, rjo) 
is called scale invariant, or a Harrizon-ZeFdovich spectrum, if the variance of the mass 
perturbation on horizon scale Ih = arj is constant, time independent: 

)/H> = COI1St - 

Here ( ) denotes statistical average over many realizations of perturbed universes with 
identical statistical properties, and we assume that it can be replaced by a spatial average. 

{5M)i H = [ 5pd 3 x = p [ (5p/p)d 3 x » pV H f " 5 k e ikx d 3 k . 
Jv H Jv h Jo 

Here we have assumed that perturbations on scales smaller than In {k > ku = I/77) 
average to zero due to the integration over Vh, and that perturbations on scales larger 
than Ih (k < ku) are approximately constant in Vh, such that integration over Vh just 
gives rise to a factor Vh- Using M = pVn, we obtain 

/•ku 

{5M/M) lH « / 5 k e lkx d 3 k , (59) 







such that 



Un[/M)T n ) * I d 3 x d 3 k £ H d 3 k'5 k 5 k ,e i( - k+k ' )x 



2j3i, _ / n/i„ „.\j3; 



\5 k \ z d 6 k = / P{k,ri)d 6 k , 
Jo 

For the first equal sign we have evaluated the spatial integration leading to a 5-distribution 
and the k! integration using 5- k = 5 k . If we now assume P(k, rf) behaves like a power law 
in k, P(k,r/) oc k m , scale invariance requires m = —3, such that 5 k (r/ k ) oc k~ 3 / 2 , where 
as before r\ k denotes the conformal time at which the scale k~ l enters the horizon, i.e. 
r\ k = k^ 1 . (We should not be concerned about the logarithmic singularity for k — > 0. 
We anyway have to introduce a low end cutoff to ( |59[ ) at about kn Q , since perturbations 
on scales larger than the present Hubble radius cannot be distinguished from background 
contributions.) 

For a scale invariant spectrum we thus have to require 

5 k {rf) oc k~ 3 l 2 on super horizon scales. 
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Since we always can choose a gauge such that Sk(rj) ~ (rjk) ^kiv)) this yields 
^k(Vk) « tffc(%) « A:" 3 / 2 . 



As is well known from linear perturbation theory (see e.g. ||24||), matter perturbations on 
scales which enter the horizon during the matter dominated era subsequently start growing 
according to Sk{f]) oc a{rj) oc rj 2 . Perturbations on smaller scales which enter the horizon 
already in the radiation dominated era, remain approximately constant until the time of 
equal matter and radiation, rj eq , and then also grow according to 5 oc a. We therefore end 
up with the following behavior of density fluctuations at late times rj > rj eq : 

, j Dk l / 2 if for Ar 1 > rj eq 

dk{rj) * \ Dk-^(r,/ Veq f for k- 1 < rj eq , 



or 



D 2 krf _ . 

(r? ' ' ^ (1 + (k/k e )2)2 ' Where ke i = r?e "' 



D is a constant of proportionality which is determined by the value of S k at horizon 
crossing, 5k(rjk) = Dk~ 3 ^ 2 . On large scale we then find 

P(i]o, k) oc k n , with n = 1. 

The power n in this law is the so called spectral index. 

We have thus seen, a power spectrum P(r]o, k) oc k n on large scales leads to the behavior 

*( Vk ,k)^k n / 2 - 2 . (60) 

Therefore, the behavior of ^>(r]k,k) in our model yields spectral indices 

_ f 4 for 5 = 

n |0 for B > (lOMpc)- 2 w (lO^eF) 2 , 1 J 



As is discussed in detail in |22|], the COBE experiment is still compatible with spectral 
indices within the bounds —2 < n < 2, with a peak in the likelyhood around n ~ 1.3. The 
case B = is therefore excluded by the COBE DMR limits on the spectral index. The 
possibility B ^ cannot be ruled out from the COBE data alone. The relatively strong 
clustering of galaxies on large scales (see e.g. even strongly disfavors a large spectral 
index, but this always depends on the matter model assumed (CDM is only consistent 
with n < 1, whereas HDM prefers larger spectral indices). 

Let us now investigate the amplitude of the induced perturbations. The value of 
e = H 2 / {-nm 2 pl ) can be related to the amplitude of the quadrupole of microwave back- 
ground anisotropies. C2 can quite generally determined by the density perturbations 
today according to (see p5[) 



c 2 ~ / \5 k (m)\ 2 {h{y)/yfdy , (62) 

Jo 

with y = kr] . Using 6 k (rj ) = ey~ 2 d (rj k) 2 / '3 = (e/3)(7] /r] end ) 2 for large enough scales 
which are relavant for the quadrupole anisotropy, we obtain in the case B ^0 
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With the COBE quadrupole of C 2 ~ 3 • 10 _n this yields H « 0.02Gey. A scalar field 
mass of m = \ B <C H , certainly incompatible with the standard model of particle physics 
which is well tested at these low energies. 

For B = one finds C2 ~ e 2 = (H/m p i) 2 , leading to H as 10 16 GeV, a typical 
GUT scale, but this possibility is already ruled out due to its large spectral index. Cer- 
tainly, as soon as reliable intermediate measurements of the cosmic microwave background 
anisotropies are available, the spectral index n will become even more constrained. 



6 Conclusions 

In this work we have investigated inflationary models in which the perturbations of second 
order in 5(p are the dominant ones. We have studied the evolution of these perturbations. 
Depending on the mass of the scalar field, we found that the resulting density fluctuations 
have spectral index n = 4 for massless fields and n = for massive scalar fields, the second 
of which is compatible with the COBE DMR results. However, the quadrupole amplitudes 
for the two cases are completely different: 

^ / (H/lGeVf for B £ 
62 ~ \ (H/m pl ) 2 for B = 0. (64) 

In the case B / 0, for an infationary model to yield small enough perturbations, one 
therefore must make sure that the perturbations of second order in the scalar field discussed 
in this work are substantially suppressed. Otherwise, due to the enormous increase cx y~ 2 d 
of the Bardeen potentials during the radiation era, they will soon dominate and spoil the 
perturbation theory. 

For 'intermediate' masses, m = \^B ~ W^^eV, the spectral index can turn out to 
be of order n ~ — 2 and the amplitude C2 ~ e 2 . (When the relevant physical scales, 
lOOMpc < k~ 1 < lOOOMpc happen to lay at the bending point of shown in Fig. 2.) But 
this requires severe fine tuning of the scalar field mass and is thus not investigated any 
further. 

We have not followed through the reheating period, the transition from the inflationary 
to the radiation dominated phase. We just assumed that perturbations induced during 
this short time can be neglected on super horizon scales, such that the Bardeen potentials 
pass continuously over the transition. If our scalar field represents the inflaton field, this 
assumption has to be justified by further work. But if we regard it as a scalar field in 
an external de Sitter universe which later becomes radiation dominated, our analysis is 
certainly valid. 

In the radiation era we have only analyzed one mode of perturbations of the scalar field. 
In principle we would have to perform the convolution presented in Section 4 to obtain 
the entire contribution to a given mode k\ of the Bardeen potentials which are quadratic 
in the scalar field. Repeating the calculation outlined in Section 4, one immediately sees 
if Ak ~ const, (which is the case for B = 0) also the Bardeen potentials turn out to be 
constant and of the same order of magnitude as for one mode. On the other hand, for 
Ak is a rapidly oscillating function of sfa. oc 1/k 2 and thus the main contribution comes 
from the wave number pair k = k\/2 and we do not expect the behavior of the spectrum 



to differ substantially from the one mode result. Therefore, we expect the results (61) and 



(63) obtained in a one mode analysis to remain valid for the full spectrum of fluctuations. 
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Figure Captions 



Fig. 1 The amplitude of the scalar field, Ak is shown as a function of y for 3 characteristic 
values of a. 

Solid line: a = 10 5 , with B/H 2 = 10" 4 , such that y end = 0.056. 
Dotted line: a = 120, with B/H 2 = 10^ 4 , such that y end = 0.03. 
Dashed line: a = 1, with B/H 2 = 10~ 4 , such that y en d = 0.1. 

Fig. 2 The amplitude of the Bardeen potential \I> at second horizon crossing, k = rf x 
is shown in arbitrary units as function of yend = 

k/(a en dH). The dashed line shows the 
k~ 2 behavior. In this example B/H 2 = 10 -4 , so that oiy 2 nd = 1 for y en d = 0.01. The 
calculated curve even becomes steeper at somewhat larger values of ay 2 and then, around 
a = 1 (corresponding to y en d = 0.1), bends over to a constant. 
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This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/astro-ph/9404043vl 



This figure "figl-2.png" is available in "png" format from: 



http://arXiv.org/ps/astro-ph/9404043vl 



